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Abstract 

A dependence of the isobaric ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of the recon- 
structed quasiprojectile for the reaction of 28 Si beam with 112 ' 124 Sn targets at two 
different projectile energies of 30 and 50 MeV/nucleon is presented. We demon- 
strate a linear dependence of the observable ln(Y( 3 H)/Y( 3 He)) on the N/Z ratio 
of the quasiprojectile and show the dependence of the slope on the reconstructed 
excitation energy of the quasiprojectile. We relate this slope dependence at a given 
excitation energy to the temperature of the fragmenting system. Using the model 
assumptions of the statistical multifragmentation model, a method of temperature 
determination is proposed. A caloric curve is constructed and compared to the result 
of the double isotope ratio method for the same set of the data and to the results 
of other studies. 
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Multifragmentation of highly excited nuclei has been subject of significant 
interest for many years. Since the nuclei are two-component systems consist- 
ing of protons and neutrons, the influence of isospin on the disassembly of 
hot nuclei is a question of general interest. Experimental studies have demon- 
strated the influence of the isospin of projectile and target on the isotopic 
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yields of fragments [1-5]. Theoretical framework for fragmentation of asym- 
metric systems have been developed [6-11] and a possible phase transition 
into an isospin-symmetric liquid and an isospin-asymmetric gas phase have 
been suggested. In our recent work [12] we presented the characteristics of the 
quasiprojectiles with Z = 12 — 15, reconstructed from the fully isotopically re- 
solved fragments with Zf < 5, detected in the forward angles using the FAUST 
multidetector array [13] at Cyclotron Institute of Texas A&M University. The 
distributions of the quasiprojectile velocity, mass, charge and excitation energy 
were obtained in the reaction of a 28 Si beam with 112,124 Sn targets at projectile 
energies 30 and 50 MeV/nucleon. We demonstrated that these experimental 
observables of the quasiprojectile can be reproduced satisfactorily by hybrid 
simulations, using a code implementing the model of deep inelastic transfer 
[14] for a calculation of the properties of excited quasiprojectiles and the SMM 
code [15] for a description of multifragmentation of excited quasiprojectiles. 
Further details can be found in [12]. The level of an agreement obtained im- 
plies that the influence of the non-statistical effects caused by violent collisions 
of the projectile and the target on the properties of detected quasiprojectiles 
can be practically excluded for the selected events. The obtained set of data 
allows further study of thermal quasiprojectile multifragmentation, especially 
the study of the influence of the quasiprojectile isospin on the properties of 
fragmenting system. This data is of specific interest because the isospin of 
the system which actually undergoes multifragmentation is known with good 
precision. The first part of the study was presented in our recent work [16] 
where multiplicities and N/Z ratios of light charged particles ( LCPs ) and 
intermediate mass fragments ( IMFs ) were studied. Different trends of multi- 
plicities and an inhomogeneous distribution of isospin between light charged 
particles and intermediate mass fragments were observed. The multiplicity of 
LCPs increases rapidly with proton excess, while the multiplicity of IMFs in- 
creases slowly with increasing neutron number. For quasiprojectiles with large 
proton excess an inhomogeneous distribution of isospin into proton-rich LCPs 
and more symmetric IMFs occurs. Furthermore, dependences of the isobaric 
ratio Y( 3 H)/Y( 3 He) on the N/Z ratio ( as a variable closely connected to 
isospin ) and on the excitation energy of the reconstructed quasiprojectiles 
were presented in [16]. In this work we will present a detailed study of the 
trends of the isobaric ratio Y( 3 H)/Y( 3 He). The relationship of Y( 3 H)/Y( 3 He) 
to thermodynamical observables will be discussed. 

The dependences of the isobaric ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of the 
quasiprojectile ( N/Zq P ) are shown in Fig. 1 for the reaction of a 28 Si beam 
with 112 ' 124 Sn targets at projectile energies of 30 and 50 MeV/nucleon. The 
quasiprojectiles with Z — 12 - 15, reconstructed from the fully isotopically 
resolved charged fragments with Zf < 5, are selected ( for details of the 
experiment see [12,17] ). Various symbols represent the experimental data for 
different targets at each beam energy. The data within each set form nearly 
straight lines thus implying an exponential dependence of the isobaric ratio 
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Y( 3 H)/Y( 3 He) on N/Zq P . For a given projectile energy, the Y( 3 H)/Y( 3 He) 
dependence on N/Zqp is practically identical for the two different targets 
ii2,i24g n rj^jg g^Q-ws that the isotopic trends are determined by the isospin of 
the excited system which fragments and the trends obtained by comparison 
of the data for different projectile-target systems describe only the overall 
trends. A detailed discussion of this issue is in our previous work [16]. For the 
two different projectile energies 30 and 50 MeV/nucleon the Y( 3 H)/Y( 3 He) 
dependence exhibits a different slope. In the above mentioned work [12], the 
distributions of the apparent excitation energy ( E* pp ) are practically identical 
for the different targets at a given projectile energy and shift to the increasing 
values of E* with an increasing beam energy. The mean values of E* app are 
101.2 and 102.3 MeV for 112 Sn and 124 Sn targets at the projectile energy 30 
MeV/nucleon and 142.8 MeV for both targets at the projectile energy 50 
MeV/nucleon. Such a similarity may imply that the slope of the dependence 
of isobaric ratio Y( 3 H)/Y( 3 He) on N/Zq P correlates with the mean excitation 
energy of the quasiprojectile. Preliminary indications of this dependence have 
been reported in ref. [16]. This assumption may be examined by extracting 
the isobaric ratio Y( 3 H) /Y( 3 He) for different bins of the apparent excitation 
energy of the quasiprojectile. 

Fig. 2 shows the dependences of the isobaric ratio Y( 3 H) /Y( 3 He) on N/Zqp for 
the nine bins of the apparent excitation energy per mass unit of the quasipro- 
jectile ( e ap p ). The data for both targets and projectile energies were com- 
bined to increase the statistics. This is possible because for a given excitation 
energy bin, the Y( 3 H)/Y( 3 He) dependences agree within the statistical devi- 
ations. The squares represent the experimental data and the lines are the fits 
obtained using the MINUIT minimization package [18]. As one can see on Fig. 
2, linearity is the overall feature of the logarithmic plots of the Y( 3 H)/Y( 3 He) 
ratio in all excitation energy bins and is especially significant in the excitation 
energy bins with a high statistics. The slopes are steepest at low excitation 
energies and become flatter with increasing excitation energy. 

Since the experimental dependences of the isobaric ratio Y( 3 H)/Y( 3 He) on 
N/Zqp correlate with e app and are linear on a logarithmic scale ( and thus 
depend exponentially on N/Zqp ), it would be of interest to relate the slope 
to the thermodynamical observables of the quasiprojectile multifragmentation. 
As a starting point we chose the statistical multifragmentation model ( SMM ) 
[15] which was employed in the simulation successfully describing our data 
[12,16]. In the macrocanonical approach, the SMM gives the mean number of 
the emitted fragments of a given N k and Z k as 



where T is the freeze-out temperature, V is the freeze-out volume, gN k z k is 



< M N k Z k >~ 9N k Z, 



• fc -^ 3/2 -p( 



F Nk z k (T, V) - jj, n N k - fi p Z k 
T 



) (1) 
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the degeneracy of the fragment ground state, Vf is the free volume available 
for translational motion of the fragment, A^ is the thermal wavelength of the 
nucleon, F NkZk (T,V) is the internal free energy of the fragment and fi n and 
fip are the chemical potentials for neutrons and protons. The assumption of 
chemical equilibrium is not incorporated into SMM [15] and the chemical po- 
tentials fi n and fip are used as a Lagrange multipliers ensuring the conservation 
of the mean neutron and proton number. When considering the isospin depen- 
dence of the experimental Y( 3 H)/Y( 3 He) ratio for a given excitation energy 
bin, the exponential dependence suggests a possibility to determine a single 
temperature for each bin as a simplest possible interpretation. The isobaric 
ratio Y( 3 H)/Y( 3 He) then can be determined as 

Y( 3 H) <M 21 > F 21 (T,V)-F 12 (T,V)-fj, n + fj,^ 

Y( 3 He) ~ < M12 > ~ 6XP T > [Z) 



where we used the fact that V) and Ay are the same for the whole partition 
and that the ground state spins of both fragments are equal. 

We approximate the fi n — fi p by the difference in the neutron ( S n ) and proton 
( S p ) ground state separation energy. According to the model of a nuclear 
Fermi gas at zero temperature, S n and S p can be identified with the chemical 
potentials for protons and neutrons in the ground state. When introducing 
finite temperature, the Fermi gas model predicts only a moderate change of 
the chemical potential for temperatures well below the Fermi energy. Thus the 
approximation of chemical potentials by separation energies can be still con- 
sidered valid for highly excited thermally equilibrated nuclei and determines 
the correct chemical potentials prior to multifragmentation. During multi- 
fragmentation the system expands which creates two different trends. When 
considering the system as Fermi gas expanding to the freeze-out density ( typ- 
ically 0.4p for the present case according to the formula used in the SMM ), 
the values of the Fermi energy decrease by more than 40 %. The same behav- 
ior can be expected also for \i n — fi p . On the other hand, chemical potentials 
of nucleons are influenced by the process of multifragmentation. As reported 
in our previous work [16], the IMFs are more isospin-symmetric than light 
charged particles. Within statistical multifragmentation model, the isobaric 
ratio Y( 3 H)/Y( 3 He) can be identified with the ratio of mean densities of free 
neutrons and protons for a given fragment partition ( as one can see after 
dividing equation 1 by freeze-out volume ). As follows from Figs. 1,2, such 
ratios differ dramatically from the starting values and one can expect also the 
values of fi n — fi p significantly higher than those of an expanded Fermi gas. 
When assuming that both effects have similar magnitude, the approximation 
used can be considered as a reasonable estimate of the actual values. For the 
full set of data used in Fig.2, fi n — fi p can be approximated as 

Hn - Hp ~ S p - S n = (-43.44 ± 0.97) + (38.57 ± 0.9A)N/Z QP [MeV]. (3) 
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Experimental mass excesses [19] have been used for the evaluation. The de- 
pendence of S n — S p on the N/Z ratio of the quasiprojectile can be considered 
linear with good precision. The difference in the internal free energy of the 
fragments 3 H and 3 He, according to the evaluation formulas given in [15], is 
of the order of a few hundred keV on a broad range of N/Zq P and changes 
much more slowly than fi n — /j, p . In this paper we neglect its weak dependence 
on N/Zqp and treat it as a constant independent on N/Zqp. Within our ap- 
proximation, the resulting expression for the isobaric ratio Y( 3 H) /Y( 3 He) will 
be 

ln(Y( 3 H)/Y( 3 He)) = \n(K(T)) + - /i p )/T (4) 



where K(T) is a proportionality factor dependent on the temperature but 
independent of the N/Z ratio of the fragmenting system. According to the 
evaluation formulas for the free energy given in [15], one may expect the 
values of proportionality factor K(T) to be close to unity. 

Since both ln(Y( 3 H)/Y( 3 He)) and \i n — \i v depend linearly on N/Zq P ( as 
may be seen from Fig. 2 and equation 3 ) and since we choose the value of 
\n(K(T)) to be a constant, independent of N/Zqp, the expressions on both 
sides of equation 4 are in fact first order polynomials of N/Zq P . Then the 
solution is obvious, constants and coefficients of first order ( slopes ) should 
be equal on both sides of the equation. This allows a determination not only 
of the temperature ( as a ratio of the slopes ) but also of \n(K(T)) from 
the comparison of the zero order coefficients ( constants ) using the extracted 
temperature. The resulting values of T and \n(K(T)) are given in Fig. 3 and 
Fig. 4 for the different quasiprojectile excitation energy bins. The numerical 
values of T and \n(K(T)) are given in the Table 1 along with the parameters 
of the linear fits of ln(Y( 3 H)/Y( 3 He)) given in Fig. 2. The dependence of the 
temperature on the excitation energy ( caloric curve ) given in Fig. 3 ( solid 
squares ) is compared to the experimental caloric curve obtained for the same 
set of the data by the double isotope ratio method [20,21] for the thermometer 
d,t/ 3 He, 4 He ( the dashed line indicates the double isotope ratio temperature 
and the solid lines indicate the statistical errors ). The formula 

Y( 2 R)/Y( 3 B.) 

T d , t/ 3 H e,,He = 14.31MeV/ ^(1.59 ^^^ ) (5) 



was used for the calculations. The agreement between the two plots is reason- 
able. Both methods give the value of the temperature between 5 and 7 MeV 
for the apparent excitation energies above 5 MeV/nucleon. Similar values of 
the temperature are reported by other experimental studies [22-25]. The val- 
ues of the temperature predicted by the SMM [15] are also in the same range. 
The values of \n(K(T)) are between and 1 and suggest the values of K(T) 
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between 1 and 3, which are in good agreement with the values of the SMM. 
The dependence of \n(K(T)) on the temperature is quite weak, which is also 
consistent with the expectations. Large relative errors of the \n(K(T)) values 
are caused by the subtraction of numbers of similar magnitude. The agree- 
ment between temperatures determined by two different methods shows that 
the assumptions made for ji n — ji p reflect the physical trends which take place 
in the freeze-out configuration. 

Since the neutron emission from the quasiprojectile is not included in the data 
presented here, we tried to estimate its influence. Using a backtracing proce- 
dure described in [12], we estimated the initial values of the excitation energy 
of a hot quasiprojectile at the stage where it also includes the neutrons that 
are later emitted. The mean value of the excitation energy taken away by a 
neutron emission ranges from 0.5 MeV/nucleon for the highest bin of e app to 
1.5 MeV/nucleon for the lowest bin. This applies to both methods of the tem- 
perature determination presented here and implies that the events represent 
a still broader range of the initial excitation energies of the quasiprojectile. 
Since the multiplicity of the charged fragments decreases with decreasing e app , 
the decrease of the determined temperatures in the lowest bins of e app may 
be a signature of the onset of a low energy deexcitation mode where several 
neutrons and light charged particles are emitted prior to the breakup of the 
partially cooled residue into two massive fragments ( such a mass distribution 
is observed in the channels with 3 and 4 charged fragments ). The temperatures 
determined thus become mean values representing the range of temperatures 
at which hydrogen and helium isotopes are emitted during the deexcitation 
cascade. The use of 3 He for the evaluation of the isobaric ratio Y( 3 H) /Y( 3 He) 
could be affected by the so-called 3 He-puzzle ( see [26] ), nevertheless it ap- 
pears that the 3 He-anomaly is not relevant for peripheral collisions induced 
by light projectile nuclei. 

We presented the dependence of the isobaric ratio Y( 3 H)/Y( 3 He) on N/Zq P 
for the reaction of a 28 Si beam on 112 ' 124 Sn targets at two different projec- 
tile energies 30 and 50 MeV/nucleon. We demonstrated a linear dependence 
of ln(Y( 3 H)/Y( 3 He)) on N/Zq P and the dependence of the slope on the ex- 
citation energy of the quasiprojectile. Using the model assumptions of the 
statistical multifragmentation model, we demonstrated the possibility of re- 
lating the slope of the dependence of ln(Y( 3 H)/Y( 3 He)) on N/Zqp to the 
temperature at a given excitation energy. We constructed a caloric curve and 
demonstrated its equivalence to the results of the double isotope ratio method. 
Our results show that the temperatures extracted from the isospin dependence 
of Y( 3 H)/Y( 3 He) are consistent with the temperature of the freeze-out config- 
uration where thermal equilibration occurs. The approximation used for the 
evaluation of the difference of neutron and proton chemical potentials is con- 
sistent with significant changes of densities of neutrons and protons during 
multifragmentation. 
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Fig. 1. Dependence of the isobaric ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of 
the fully isotopically resolved quasiprojectiles with Z tot = 12 — 15. Squares - 
28 Si(30MeV/nucleon)+ 112 Sn, circles - 28 Si(30MeV/nucleon)+ 124 Sn, up triangles - 
28 Si(50MeV/nucleon)+ 112 Sn, down triangles - 28 Si(50MeV/nucleon)+ 124 Sn. 
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Fig. 2. Dependence of the isobaric ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of the 
fully isotopically resolved quasiprojectiles with Z to t = 12 — 15 for the different bins 
of € app . The panels (a) - (i) correspond to the excitation energy bins given in Table 
1. The full set of the data from reactions of 28 Si beam with 112,124 Sn targets at the 
projectile energies 30 and 50 MeV/nucleon was used. 
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Fig. 3. Dependence of the temperature T, determined from the dependence of the 
isobaric ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of the fragmenting quasiprojectile, on 
e app ( solid squares ). Dashed and solid lines indicate the values and the statistical 
deviations of the temperature determined using the double isotope ratio method for 
the thermometer d,t/ 3 He, 4 He. 
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Fig. 4. Dependence of ln(K(T)), determined from the dependence of the isobaric 
ratio Y( 3 H)/Y( 3 He) on N/Z QP , on e app . 
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Table 1 

Parameters of the exponential fit y(3}{e) = ex P( ai + a 2("f")<2-P) f° r the different bins 
of the apparent excitation energy e app and the determined values of the temperature 
T and Yn(K(T)) for the same bins. 



€ ap p, MeV/nucleon 


ai 


02 


T, MeV 


ln(K(T)) 


2.5-3.5 


-17.895±2.170 


15.719±2.007 


2.45±0.43 


-0.19±3.19 


3.5-4.5 


-12.665±1.175 


11.642±1.160 


3.31±0.44 


0.45±1.81 


4.5-5.5 


-9.124±0.517 


8.730±0.538 


4.42±0.38 


0.71±0.86 


5.5-6.5 


-7.267±0.483 


7.099±0.515 


5.43±0.53 


0.73±0.80 


6.5-7.5 


-6.808±0.518 


7.031±0.564 


5.49±0.58 


1.11±0.86 


7.5-8.5 


-5.634±0.682 


5.868±0.744 


6.57±1.08 


0.97±1.10 


8.5-9.5 


-5.023±1.084 


5.218±1.210 


7.39±2.16 


0.85±1.75 


9.5-10.5 


-6.499±2.020 


6.846±2.202 


5.63±2.30 


1.21±3.21 


10.5-12.5 


-5.323±2.931 


5.689±3.241 


6.80±4.86 


1.08±4.69 
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